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Abstract 

^^ . It is demonstrated that a nonrelativistic quantum scale anomaly manifests itself in the appear- 

_h \ ance of composite operators with complex scaling dimensions. In particular, we study nonrelativis- 

tic quantum mechanics with an inverse square potential and consider a composite s-wave operator 
O = iptp. We analytically compute the scaling dimension of this operator and determine the prop- 
agator (0\TOO'\0). The operator O represents an infinite tower of bound states with a geometric 



Q\ 



Q_i' energy spectrum. Operators with higher angular momenta are briefly discussed. 
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I. INTRODUCTION 

Relativistic conformal field theory is a well-developed subject with numerous applications 
in statistical physics (systems near a continuous phase transition) and high-energy physics. 
Recently, Nishida and Son made an important step towards systematic understanding of 
nonrelativistic conformal field theories (NRCFT) [l|, I2J 1 defined as invariant with respect to 



Q. 



the Schrodinger symmetry [4j. This symmetry contains the usual Galilei symmetry of nonrel- 
ativistic systems extended by a scale and special conformal transformation. The symmetry 
transformations form a group, called the Schrodinger group, which is a direct nonrelativistic 
analogue of the relativistic conformal group. As a simple example, the free nonrelativistic 
field theory has the Schrodinger symmetry. However, also a number of theories with inter- 
actions such as cold fermions at unitarity are believed to be invariant with respect to the 
Schrodinger symmetry, which provides powerful constraints on the correlation functions. 
Similar to relativistic conformal field theories the basis of a NRCFT is formed by primary 



operators 



2J. A local primary operator 0(t,x) has a well-defined scaling dimension A^ 



and particle number Nq 

[D,0} = iA 0, [N,0} = N 0, (1) 

where O = 0(t = 0,x = 0), D and N denote the generators of scale and particle num- 
ber symmetry, respectively. In addition, the primary operator O must commute with the 
generators of Galilei boost K t and the special conformal generator C 

[K u O] = 0, [C,O] = 0. (2) 

The set of descendant operators constructed by subsequent application of time and spatial 
derivatives on the primary operator O form the irreducible representation of the Schrodinger 
algebra |l|. The scaling dimensions of descendants are related in a simple way to the scaling 
dimension of the parent primary operator O. 

It is known that in some nonrelativistic theories the classical Schrodinger symmetry can 
be broken by a quantum scale anomaly 5j. 2 This often leads to discrete scale invariance and 



1 for an earlier work see 

2 In general, a quantum anomaly means that a classical symmetry is broken at the quantum level due to a 
regularization and renormalization procedure. 



a geometric bound state energy spectrum 6|, |7|] (for a review and source of further references 
see 8j). The anomaly also manifests itself in the renormalization group flow of a contact 
coupling which develops a limit cycle [8j. Another signature of the nonrelativistic scale 
anomaly is the appearance of composite operators with complex scaling dimensions. In this 
paper we concentrate on the latter and discuss some properties of composite operators with 
complex scaling dimensions in a nonrelativistic quantum theory. In particular, we consider 
a two-particle problem in the quantum theory with an inverse square potential interaction. 
In this simple example we analytically compute the propagator of the s-wave composite 
operator O = ipip and examine its structure. 

The rest of the paper is organized as follows: In Sec. [IT] we provide a brief introduction to 
quantum mechanics with an inverse square potential. In Sects. IHIirVl we consider the s-wave 
composite operator O = ipip and calculate its two-point correlation function in momentum 
space. Sec. [Vljis devoted to a short discussion of composite operators with higher angular 
momenta. Finally, we close in Sec. IVHI with conclusion and discussion. Technical details 
are summarized in Appendices. 

II. QUANTUM MECHANICS WITH INVERSE SQUARE POTENTIAL 

Consider the nonrelativistic quantum theory of identical bose particles interacting 
through a long-range inverse square potential V(r) = — -% in D spatial dimensions. 3 The 
corresponding nonrelativistic quantum field theory is defined by the microscopic (classical) 
action 

S[^^*]= f dt f d D xip*(t,x)[id t + —}ip(t,x) 

Iff 2 (3) 

--J dt J d D xd D y r(t,x)r(t,y)v(\x-y\Mt,yMt,x). 

In our convention the particle mass M^ and the reduced Planck constant H are set to unity. 
The classical action fl3]) is symmetric under the global U(l) internal transformation, associ- 
ated with the particle number conservation, and is invariant under the Galilean spacetime 
symmetry group. The interaction parameter k characterizes the strength of the long-range 



3 In this work we use field-theoretic (second-quantized) formulation of nonrelativistic quantum theory, which 
allows to treat different few-body problems in a unified way. 



inverse square potential and is positive in the attractive case. Most remarkably, in any 
spatial dimension the action ([3]) is scale invariant because the potential is a homogeneous 
function of degree —2 and has the same scaling as the kinetic part in Eq. (EJ). 

The two-particle sector of quantum theory with the inverse square potential was studied 



extensively |5|, I9M12J and by now is well understood. It is a paradigm for nonrelativistic 

conformal invariance and scale anomaly with a number of physical applications ranging 

rom low-energy atomic to high-energy particle physics (for the list of applications see e.g. 



12|). 

The Feynman rules for the particle propagator iG^ = (0|T^M|0) and the interaction 
vertex in momentum space are 



p 2 _i_ • 




= 2iKF D (l\ 



(4) 

where ie is an infinitesimal imaginary number, ensuring the retarded causal structure of 
the propagator. I = p 2 — p i = p^ — p^ is the spatial momentum transfer during a collision 
(/ = |/|), and Fd(1) denotes the Fourier transformation of the inverse square potential which 
in D > 2 reads 

F D{ 1) = / A± e*p H f. A = W"W-l)P-° , (5) 

In the rest of this work we restrict our attention to spatial dimension D > 2. 

III. COMPOSITE OPERATOR O = iptp 

Consider a local two-particle s-wave operator 0(t,x) = ipip(t,x) which annihilates two 
identical bose particles at the spacetime point (t,x). As demonstrated in Appendix |Aj we 
can assign to the operator O a pair of scaling dimensions 



D + 2 (D-2) 2 

The value of the interaction coupling k determines two qualitatively different regimes: 



2 



Undercritical regime k < n cr = - ~ — : Both A + and A_ are real and the action 
((3]) defines two conformal quantum field theories CFT + and CFT_. More precisely, 
the quantum field theory defined by the classical action (J3J) flows towards CFT + in 
the infrared and towards CFT_ in the ultraviolet in the renormalization group sense. 
To obtain the action of CFT + (CFTJ), we must add to the inverse square potential 
a delta function localized near the origin 

V{r)^~--8{r-e) (7) 

r A e 



and fine-tune the value of the dimensionless contact coupling A [13] to the infrared 
(ultraviolet) fixed point. The operator O is a nonrelativistic two-particle primary 
operator because it is a local product of two primary operators ip. 

Overcritical regime k > n cr : The quantum field theory ceases to be conformal and 
renormalization group evolution of the coupling A in Eq. (0 develops a limit cycle 



|allfiUl2 



12|. A± in Eq. (jHJ) becomes complex and form a conjugate pair. Due to the 
loss of conformality A± do not fulfill the definition for scaling dimension as given 
in Eq. (JT]). However, we will still use the term scaling dimensions for A± even in 



the overcritical regime. This is motivated by observation made in 12] that one can 
formally extend parameter space of the contact coupling A to the complex plane. In 
this extension two complex fixed points emerge in the overcritical regime. At these 
fixed points the operator O has scaling dimensions A± in the sense of Eq. ([T]). 

We now calculate the two-point function (two-particle propagator) 

iG (t 2 ,£2]ti,xx) = (OlTC^^O^xOlO) (8) 

in the overcritical regime k > K cr . Intuitively, the two-point function is proportional to the 
probability of creating two identical particles at the spacetime point (£i,xi) and subsequent 
destruction of the pair at the distinct point (£2, #2)- 

The translational invariance of the action ([3]) allows us to transform %Gq to the momen- 
tum space 

iG (u, p) = f dtd D xe iuJt - l ^iG (£, x\ 0, 0) . (9) 

The diagramatic expression of the propagator zGo(a;,p) is depicted in Fig. [U The first 
diagram in Fig. [1] corresponds to a free propagation of two particles, while the remaining 
sum of diagrams takes into account the influence of the long-range interaction. 



iG (w,p) = 



Figure 1: The two-particle propagator iGo(tJ,p) illustrated as a sum of Feynman diagrams. The 
full circles denote the composite operator insertions. 

At this point it must be stressed that few-body problems in nonrelativistic quantum field 
theory are relatively simple because only particles (but no antiparticles) exist as excitations 
of the nonrelativistic vacuum. 4 This remarkable fact and the particle number conservation 
lead to the important diagramatic simplification: Any loop with line arrows pointing in the 
same direction vanishes. Mathematically, this can be demonstrated with the help of the 
residue theorem for the frequency loop integration. For this reason the above-mentioned 
loops which involve both particles and antiparticles vanish. 

The diagramatic simplification has a number of important consequences: 

• The particle propagator is not renormalized and coincides with the bare propagator 
depicted in Eq. (J3J). For this reason there are only bare propagators in Fig. [TJ 

• The long-range interaction vertex is not renormalized (i.e. there is no screening). That 
is why there are only bare vertices in Fig. [TJ 

• There are no "crossed" interaction lines in Fig. [TJ 

We note that due to the Galilean symmetry of the action ([3]), the two-particle propagator 
iGo(u,p) should be a function of only the Galilean invariant combination u — p e /4. Hence, 
it is advantageous to switch to the center-of-mass frame (i.e. set p — 0) first and calculate 
iGo(u) = iGo(u,p = 0). In the very end we can recover momentum dependence by the 
substitution iGo(u) — i iGo(u — p 1 /^). 

Separation of relative and center-of-mass motion in the two-body problem allows us to 
relate the two-particle propagator iGo{oj) to the energy Green function Gd{t" \f';ui) of a 



4 From this perspective, many-body problems are more difficult due to a more complicated vacuum state 
which allows the presence of hole excitations. 
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iGoM = j r „ k {— £-*('"* - ?) + e~ i nr=i /r ( _tfc£t} 

= iG D (f" = 0,f' = 0;w) 

Figure 2: The relation between the two-particle propagator iGo(ui) and the energy Green function 
Gd{t" ' , r';uj). Full lines with an arrow denote here the bare propagator of a particle of the reduced 
mass M re d = i. 

single particle with the reduced mass M re ^ = —^- = | in the external inverse square potential. 
To clarify this point, consider the first diagram in Fig. [TJ which we shall call iG o (uj) 

dlo i i f f i 

TJ 2vr / - | + ieu - l - f + ie JiJv w - I 2 + ie 



j(f " = 0|(w - p 2 + ze)" 1 ^ ' = 0) = iG° D (f" = 0, f ' = 0; w) 



where Jj-= J i^ ■ In the first line we performed frequency loop integration. In the last line 
the definition of the energy Green function G° D (r", r '; a;) evaluated at the origin r " = r ' = 
was recognized. In the same way, we can perform frequency loop integrations in the diagrams 
with interaction vertices in Fig. [1] and find that 

iG {uj) = i{f" = 0|(w - H + ie)- l \r' = 0) = iG D (f" = 0,r' = 0; w), (11) 

where H = p 2 — -^. In terms of Feynman diagrams the last relation can be expressed as 
depicted in Fig. |2j The energy Green function Gd{t" ', r '; a;) for the inverse square potential 



problem was calculated in 14J, and a modified derivation is presented in the next section. 



IV. ENERGY GREEN FUNCTION G D {r" ,r';uj) 

The energy Green function in position representation in D spatial dimensions is given by 

G D (f",r';cu) = (f"\(cu - H + ier 1 \r'), (12) 

where H = p 2 + V(r) is a single particle Hamiltonian. The definition (1121) leads to the 
following inhomogeneous differential equation 

[V|„ + cu - V(f ")] G D (r", r '■ u>) = 5(f" - f '). (13) 



For a central potential Eq. (fT3|) can be solved by separation of variables. In particular, the 
energy Green function can be expanded in D-dimensional partial waves 

oo di 

G D (r",r';u;) = (rV)-^ E £ >Wn '0*^(0 ^V';"), (14) 

l=Q m=l 

where Yi m (il) denote D-dimensional spherical harmonics and d\ = ndj^v. ~ Substi- 
tution of Eq. (fl4|) in Eq. (1131) gives the differential equation for the radial energy Green 
function 



d 2 



dr 



n2 



+ u-V{r") 



{l + v?-\ 



. //2 



Gi{r",r';u)=S{r" -r') 



(15) 



with v = ^2^- Solution of this equation is now a textbook problem (see e.g. 15|). The 
radial energy Green function can be expressed as 



G,(rV';w) 



.(<)/ 



,(>)/ 



uj '{r < )uj '(r>) 
M^{«, (<) ,«, (>) }(r') 



(16) 



with r< (resp. r>) denoting the smaller (resp. larger) one of r' and r". W{u\ \u\ } 
represents the Wronskian of uf and uf, i.e. 

W{«W ) «W} = «{<M >)/ -«, (<)/ «i (>> . (17) 

The function u\ (resp. u\ ) solves the homogeneous differential equation 

uiir) = 



d 2 | u (/ + v) 

dr 2 r 2 



2 1 
4 



(18) 



with a regular boundary condition at r = (resp. r = oo). 

Consider now the inverse square potential V(r) = — \. Depending on the value of the 
coupling strength k, the solution of Eq. f fl8|) has two qualitatively different branches. We 
analyze these two regimes separately. 

A. Undercritical regime k < (I + v) 2 



Two linearly independent solutions of Eq. ( TTBl are 



Ul (r) = {yfrI Sl (£r),y/rK Sl {(r)} 



(19) 



5 We use the notation from [14| for the labelling of spherical harmonics in D dimensions. 



with C = y/—uj — ie and si = y/(l + v) 2 — k. I Sl and K Sl denote modified Bessel functions. 
We find 



«, (<) (r) ~ v^^(Cr), n; <} (r)^0. 
Substitution of Eq. ( I20l) in Eq. ( [161) yields the radial Green function 



(20) 



G l {r",r' ] u) = -V^I Sl ((r<)K Sl ({ r> ). (21) 

B. Overcritical regime k > (I + u) 2 

In this case two linearly independent solutions of Eq. ([TBI) are found to be 

m{r) = {VFl Wl (Cr), v^A^(Cr)} (22) 

with 6* = y^n — (I + v) 2 . It is straightforward to determine u\ , which is 

u { >\r) ~ VfK^iCr), ^ (>) (r) ™ 0. (23) 

On the other hand, determination of u| (r) turns out to be more subtle. It is easily 
demonstrated that any linear combination of two independent solutions ( )22l) approaches 
zero as r — > 0, and thus satisfies the regular boundary condition. Mathematically, the 
problem originates from the singular behavior of the inverse square potential at r = 0. In 
the overcritica l reg ime the potential is too singular and must be regularized. Regularization, 
introduced in 14|, imposes the boundary condition at some small but finite r = a. Here, 
following ll| . we use a different but equivalent regularization procedure which consists in 
fixing a phase angle of u\ (r) as r — > 0. Specifically, we take the function u\ (r) to be 

«j (<) ir) ~ v^ [e iri{0 h 6l ((r) + e^/^ (Cr)] (24) 

with 77(C) = — 61 In (^r)- Here we introduced a momentum scale Q that determines the 
phase shift of u\ (r) near the origin 



u) r ~ 



v^cos(^lnr + /3 + ^lnC), (25) 



where (5 = arg ( jwr^j \ ) • We emphasize that the phase shift 5i<p is a physical parameter 
that fixes u\ (r) in the unique way (up to normalization). Substitution of Eqs. (123] 125]) in 
Eq. ( 1T6|) yields the radial energy Green function 



n ( „ , x v^V 7 [e^)/ ( Cr< ) + e -^(0/ ( Cr< )] ^,(Cr>) 
Gj(r,r;w) = ;T~ T7Z i777777^ ■ ( 26 ) 



2cos(^ln(C/CD) 
V. TWO-PARTICLE PROPAGATOR iG 

We are now in the position to finish the calculation of the two-particle propagator iGo{&) 
in the overcritical (anamalous) regime. To this end, according to Eq. (II ip one must evaluate 
the energy Green function Gr>(f", r'\ u) at r" = f ' — 0. Inspection of Eqs. (I14II26P reveals, 
however, that for D > 2 this leads to a divergent result. This is not a surprise because the 
operator O is a local product of two elementary bose operators ip. Composite operators 
usually lead to additional divergences in quantum field theory which must be renormalized 
separately. We will deal with this problem in Appendices [B] and The propagator of the 
bare operator O reads 

//70 ' r /70 " 
JdZi/ ^ G °(r",?>;u)l ii=ri=r , (27) 

where S D ~ l = 2 r T p/2 x . In the last expression we introduced two angular integration aver- 
ages that should be performed before the radial limit. This ensures the final result to be 
independent of the directions of r ' and f". We now substitute the partial wave expansion 
(EH} in Eq. (ETJ) and obtain 

«Gn w =-i lmr t-— ; t- 7- NN — d , 28 

1 J »^o 2cos(^ln(C/Co*)) 

where only the s-wave (/ = 0) survives the angular integrations. Finally, we perform the 

limit r — > and find 

iG r S n (oo) = i tan U In V ~^~ ^ , (29) 

where computational details and the explicit definition of the renormalized two-particle 
propagator iG^ n (u) can be found in Appendix [Bl For the sake of completeness we compute 
iG r Q n (u) in the undercritical (conformal) regime in Appendix O 

Our final result (129]) is remarkably simple and has the following properties: 
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For negative energies u < the two-particle propagator iG r Q n (u) has an infinite num- 
ber of simple pole divergences at 

u n = -Co 2 exp (-— + -), neZ. (30) 

Hence, the composite field O* represents infinitely many different stable particles (two- 
particle s-wave bound states). The energy spectrum has an accumulation point at 
u — > as n — y oo. It exhibits geometric behavior 

— = exp \-%) ■ (31) 

For positive energies u > the two-point correlation function iG r ^ n (u) has a branch 
cut with a branch point at u = 0. 

The regularization momentum scale Q introduced in Sec. IIVI is a physical parameter 
of the theory. It determines the value of the reference (for example n = 0) two-particle 
bound state energy 

Wo = -Co 2 exp(^). (32) 

The propagator iG T Q n (u) is invariant under the discrete scaling symmetry 

Q^Qexp(j^Y meZ. (33) 

This transformation maps ui n — > uj n+m , but does not alter the energy spectrum. Phys- 
ically inequivalent Q may thus be chosen to lay within the interval ( * G [1, exp ( f- ))■ 
This supports our observation in Sec. IIVI that Q determines the phase shift angle 
#00 = #o hi Co °f the wave function near the origin. We conclude that physically 
inequivalent phase shifts <5o0 span the interval [0, 7r). 

Finally, as advocated in Sec. IHH the Galilean symmetry allows us to recover the 
momentum dependence of the two-particle propagator 



V 2 



iG r S n (oo, p) = % tan | O In -* — | . (34) 

The position space representation of the propagator is derived in Appendix [Dl 
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VI. COMPOSITE OPERATORS 00) WITH HIGHER ANGULAR MOMENTUM 

The energy Green function (TT2l) contains also information about the two-particle com- 
posite operators O^ which carry the higher angular momentum I. In the conformal regime, 
these operators are nonrelativistic primaries of ref . |l| . The p-wave and d-wave operators 
in a theory with a one-component bose field ip are explicitly given by 

Of =1) = V^-^V^, (35) 

Of= 2) = [^ViVj-V - VrfVtf] - % \4>M - V k ^V k ^} (36) 



Up to normalization these operators are fixed by two requirements. First, they must be 
constructed from only two elementary fields ip and spatial gradients. Second, they are 
nonrelativistic primaries of ref. |l| and satisfy Eqs. (0Q |2]). We expect that in a similar 
fashion the two-particle primaries with I > 2 can be constructed. It is clear from Eq. (j55|) 
that in the bosonic theory (J2J) the p-wave operator 0\ ~ vanishes, as well as all primary 
operators with odd angular momentum. The d-wave operator 0\~ is a symmetric, traceless 
tensor. 6 

The scaling dimensions of the composites 0^> were computed in Appendix |A] and are 
given by 

4? = £±i ± J(l + °^.)'-*. (37) 



They become complex for k > (l + ^-^) ■ 

The two-particle propagator iG (i) = (0\TO^O"^\0) can be extracted from the I th par- 
tial wave of the energy Green function, i.e. it is encoded in the radial Green function 
Gi{r" ,r';uj). Following closely the steps from Sec. |V] in the overcritical regime we obtain 
for the renormalized propagator 



^T~ 



iGom (w,p) = i tan [ 0, In -* | T ix „ Auh „ 4l (p), (38) 



In order to make the tensor traceless, one subtracts in Eq. (|36|) the s-wave part O s — -0A-0 — VfcV'VfeV'- 
It is instructive to express this operator as O s = — iVfcVfeO — 2idtO + 4-0 (idt + y) ip, where the first 
and the second terms are the descendants of the primary O and the third term is a primary operator. 
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where Ti 1 _,_i l j 1 _„j l is a tensor symmetric and traceless in the indices i\. . .i\ and j% ■ ■ -ji- This 
tensor does not depend on the energy u, and thus we can extract the energies of the two- 
particle bound states with angular momentum / from Eq. (13"8j) . Similar to the s-wave bound 
states, for even angular momenta I the energy spectrum exhibits geometric behavior 

u n+1 ( 2tx 



VII. CONCLUSIONS 



«. ^\-t ] - (39) 



In this work we examined composite operators in a nonrelativistic quantum field theory. 
In general, the presence of such operators is a signature of the quantum scale anomaly. As a 
concrete example, we considered quantum mechanics with the classically scale invariant in- 
verse square potential and studied the two-particle local scalar composite operator O = ipip. 
We determined the complex scaling dimension of this operator and analytically calculated 
the two-particle propagator (0\TOO*\0). In the nonconformal (anomalous) regime the op- 
erator O represents an infinite tower of s-wave two-particle bound states which form a 
geometric energy spectrum. 

For simplicity, in this work we considered the specific two-particle inverse square potential 
problem, which is the paradigmatic example of a nonrelativistic theory with the quantum 
scale anomaly. Nevertheless, we expect that our simple result f )34|) for the two-point corre- 
lator (0|TOO'|0) is universal and applies to other more complicated nonrelativistic theories 
containing composite operators with complex scaling dimensions. One prominent example is 
a system of three identical bosons interacting through a contact potential tuned to the uni- 
tarity point. Due to the Efimov effect [6|, the local atom-dimer composite operator O = ijxp 
acquires the complex scaling dimensions 

A± = ^±i So , (40) 

where the Efimov parameter s ~ 1.00624. A similar system that exhibits the quantum scale 
anomaly in the three-particle sector is the two-component fermionic system with unequal 
masses [6J]. Another interesting example is the one-dimensional model with a four-body 



resonant interaction 
is complex, too. 



16] where the scaling dimension of the five-boson composite operator 
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Recently we used the method of the nonrelativistic AdS/CFT correspondence to study 
operators with complex scaling dimensions J17J. Using the standard AdS/CFT technique, we 
calculated the two-point correlation function (0|OO'|0). The result has the same functional 
form as found in this work in Eq. (j54"|) . However, while Eq. 034|) depends on the physical 
momentum parameter Q, the holographic result of (l7| depends on the value of the UV 
momentum cutoff. The latter is unphysical, and we expect that the cutoff dependence in 
the AdS/CFT calculation can be eliminated by inclusion of a proper boundary counterterms 
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Appendix A: Scaling dimension of composite operator O = tpip 

An elegant way to calculate the scaling dimensions of the composite two-particle oper- 
ator O = ipip (and more generally of the two-particle primary operators O^ carrying the 
angular momentum I) is to employ the operator /state correspondence [l|, |l9|, |20|. To this 



end one considers two particles interacting through the inverse square potential confined 
in a harmonic trap. The total Hamiltonian H of this system can be separated into the 
center-of-mass and relative parts 



H = H R + H r , 




h r = -\v\ + Jr\ 


E° R 


^2 K u 2 r 2 
H — Y 7 1 




Ur X r r 2 1 4 ' 





D 



(Al) 



/ + ^- 2)X 



where R and r are the center-of-mass and relative coordinates. In addition, E R denotes the 
ground state energy of Hr, and E rl stands for the lowest energy of the Hamiltonian H r 
in the subspace of states with the angular momentum I. 7 According to the operator/state 



7 The relative Hamiltonian H r defines a quantum mechanical problem of a particle in the combined inverse 
square and harmonic potential, also known as the Calogero problem. The energy spectrum is formally 
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correspondence, the scaling dimension of the composite primary 0^ l \ carrying the angular 



momentum I, coincides with the lowest energy Ef of the total Hamiltonian H, expressed in 



the units of the trapping frequency lo 



(0 _Ef _D 




A^ = — ± ^ i+ — ) -«. (A2) 

Note that the operator O"' composed of two identical bose (fermi) fields vanishes if the 
angular momentum / is odd (even). This is due to the fact that two identical bosons 
(fermions) can not be in the quantum state with odd (even) angular momentum. For the 
s-wave (I = 0) operator O" =0 > = ipip we thus obtain 

A ± ^r°>^ ± ^) 2 -, (A3, 

We must stress that the operator/state correspondence applies only to nonrelativistic 
conformal field theories, and thus our result ( 1A2I) holds for k < n cr — (I + ^y^) , when both 
A^ are real. Nevertheless, it turns out that even in the anomalous (nonconformal) regime 
for k > K cr Eq. (|A2[) leads to the correct scaling dimensions A± . We illustrate this fact on 
the example of the s-wave composite operator 0(t,x). Following the observation, made in 
[I], the proper definition of the composite 0(t, x) for any value of k is given by 

0(t, x) = lim \x - y\^ip(t, x)i/j(t, y), (A4) 

where 7 is a leading near-origin power law exponent of the zero-energy wave function of the 
relative Hamiltonian H r . It can be determined from the equation Hj.r' 1 = and reads 




(-2-) 

The prefactor \x — y | 7 in Eq. (IA4|) is needed to make matrix elements of the operator 0(t, x) 
between any two states in the Hilbert space finite. From Eq. (1A4|) we can read off the scaling 
dimension of the operator O by a simple counting 

A ± = 2A^ + 7 =^±J(^) 2 -« (A6) 

which is in agreement with Eq. ( 1A3I) . found from the operator/state correspondence. 



given by two equidistant towers built on top of two "lowest state energies" E^ t and E~ t 21 1 . Physically, 
there are two proper choices of the near-origin boundary condition of the wave function that distinguish 
between + and — branches of the energy spectrum. The choice E^ { (E~j) corresponds to CFT + (CFT_) 
of Sec. [OH 
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Appendix B: Details of calculation of iG T Q n (u)) 

First, the limit r — > in Eq. (1281) can be simplified by exploiting the near-origin asymp- 
totics of the modified Bessel functions 






/., o ( x )^ fee ^olnx +/3 ) 5 fee ^ 



r(i + ^ ) 

K Wo (x) ^ a(e i{0olnx+a) + e-^oinx+a)^ ae ia = 2 - 1 -^T(-i9 ) 
with the result 



(Bl) 



D cos(#o m Co r + P) cos(6 l o In (r + a) 



iGn(oj) = — 2iab\imr 

1 J '-*> cos(0 o ln(C/Co*)) 

=2iab lim r 2 " D cos(# In Qr + P) sin(0 o In Qr + a) x (B2) 

tan ( O In — J - cot(0 o hi Co r + «) 

In the second line we used the identities ln(r = ln^r + ln((/Q) and cos(v4 + B) = 
cos A cos 5 — sin A sin B. At this point we observe that the overal normalization factor 
and the second term in the bracket of Eq. (IB2p are actually ill-defined in the limit r — > 0. 
Physically, this problem originates from the renormalization group limit cycle scaling in the 
overcritical regime. The limit r — > is intricate and will be performed in two separate 
steps. First, we subtract the energy-independent second term in the bracket of Eq. (1B2J) . 
This choice effectively prescribes the initial position on the limit cycle in the UV and thus 
fixes the energy spectrum (see Eq. f)30p ). In the second step we perform a multiplicative 
renormalization of the composite operator O by introducing its renormalized version 

Qren = j^ r V j cog ^ o j n £* f + ^ gin ^ o ln £* r + a j]-l/2 Q ^ 

with J\f = (2ab)~ . Here we emphasize that the renormalization function is not a simple 
power of r because in the overcritical regime the theory is not conformal but undergoes 
the limit cycle RG flow. After these two renormalization steps we obtain the renormalized 
two-particle propagator iG T Q n {uS) 



iG r ™(uj) = (0|TC ren O tren |0) = ztan (fl hi - " — J , (B4) 

where we substituted ( = y/—u — ie. Finally, we note that although the renormalization pro- 
cedure might appear ad hoc at first sight, in general it parallels its undercritical (conformal) 
counterpart which is described in Appendix 
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Appendix C: Two-particle propagator in undercritical regime 

In the undercritical (conformal) regime we follow directly the steps from Sec|V]and obtain 

r Act i r rjci" 
iGo(u) = (0|TOOt| ) =zlim o y A_ J J^G^rV'ju;) 






(CI) 



t\\mr 2 - u I S0 ((r)K S0 ((r), 



MO 

where in the second line we used Eqs. (lT4 t l2~lT) . Now by employing the near-origin asymptotic 
behavior of the Bessel functions fIBll) for #0 = — is and x = (r, we arrive at 

iG (uj) = -iabe i{a+p) \irar 2 ~ D ((Cr) 2s ° + e~ 2ia ) . (C2) 

The bare two-particle propagator is divergent for D > 2 and must be renormalized. This 
can be achieved in two steps, and we present the procedure for the example of CFT + . In 
the first step we neglect the energy-independent second term in the bracket of Eq. (IC2I) . 
Physically, this corresponds to probing the large distance infrared physics which is governed 
by CFT + (see Sec. IIIII) . In the second step we perform a multiplicative renormalization by 
introducing a renormalized operator 

This is necessary due to the composite nature of the operator O. Here for simplicity we 
absorbed the finite constant J\f = (— abe l ^ a+l3 '^ into the definition of the renormalized 
operator. With this definition we obtain 

iG r S n {oj) = (0\TO ren O jren \0) = i( 2so =i(-u- ie) S0 . (C4) 

Finally, by recovering the momentum dependence, we obtain a cutoff-independent, renor- 
malized two-particle propagator 

iGF(u,fi=i(-cj + £-%€} °. (C5) 

It is straightforward to transform to the position space 

J ( } , ^ (C6) 



C D 9(t)t-*+ex V [-iNo — 
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with Cd = f^zr ] (~) > ^+ = ^2^ + s o an d ^o = — 2. This is a familiar result for 
the retarded propagator of an operator O with the scaling dimension A + and the particle 
number Nq = —2 in a nonrelativistic conformal field theory. Its functional form (up to the 
normalization constant Cd) is fixed by the Schrodinger symmetry [lj]. 



Appendix D: Two-particle propagator in position space 

In this appendix we attempt to compute the two-particle propagator in the position space. 
This can be achieved via the inverse Fourier transformation of Eq. (j3] 



p- 



2vr (2n) D 



iG r S n {t,x) = / —-—^e-^+^itan # hi^ -f . (Dl) 



First we perform the angular integration and obtain 



iG r en (t, x) = ( — ) J dpp (J) J D/2 -i(px) J ^e-^z tan £ In V _i 

(D2) 



:> 



Now we introduce a new variable IV = — cu + v and get 



1 

D/2 



iG r S n {t,x) = I — J / dpp(^-J J D /2-i(px)e '*' / — e 4Vyf 2tan(^ ln J 

°(D3) 
Finally, the momentum integral can be done analytically with the result 

iG r en (t,x) = \-±) exp i-iNo^j S(t), (D4) 

where N a = -2 and S(t) = f ^e iWt i tan (e In ^f^) . We were not able to perform 
the integral over W explicitly. Based on the dimensional argument S(t) = t~ 1 f((Q) 2 t,6o), 
where / is some function of the dimensionless arguments (Co) 2 ^ an d #o- We checked that 
the function S(t) is not restricted by the Galilean symmetry. 
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